INTERPOLATION MAPS AND CONGRUENCE DOMAINS FOR 

WAVELET SETS 
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Abstract. It is proven that if an interpolation map between two wavelet sets 
preserves the union of the sets, then the pair must be an interpolation pair. 
We also construct an example of a pair of wavelet sets for which the congruence 
domains of the associated interpolation map and its inverse are equal, and yet the 
pair is not an interpolation pair. The first result solves affirmatively a problem 
that the second author had posed several years ago, and the second result solves 
an intriguing problem of D. Han. The key to this counterexample is a special 
technical lemma on constructing wavelet sets. Several other applications of this 
result are also given. In addition, some problems are posed. We also take the 
opportunity to give some general exposition on wavelet sets and operator-theoretic 
interpolation of wavelets. 

Dedicated to Larry Baggett for his great friendship, his love of 
mathematics, and his continued support of young mathematicians. 

1. Introduction 

An orthonormal wavelet is a single function in /^^(M") whose translates by 
all members of a full-rank lattice followed by dilates by all integral powers of a 
real expansive matrix on M" generates an orthonormal basis for /^^(M"). By the 
term wavelet set we mean a measurable subset E C M" with the property that the 
inverse Fourier Transform of its normalized characteristic function, !F~^( — ^yr), 
is an orthonormal wavelet. In [DL] an operator-theoretic technique for working 
with certain problems concerning wavelets was introduced that was called operator- 
theoretic interpolation. If V and r] are orthonormal wavelets in the same space, 
and ipn,! and rjn^i are the corresponding wavelet bases, then the unitary operator 
determined by the mapping ipn^i to rjn^i was called the interpolation unitary between 
tp and 7]. These interpolation operators associated with ordered pairs of wavelets 
play an essential role in the theory. They are associated with the von Neumann 
subalgebras of the so-called local commutant space, whose unitary groups provide 
natural parameterizations of certain families of wavelets. In the special case where 
the interpolation operator is involutive (i.e. has square I) the pair of wavelets is 
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called an interpolation pair of wavelets. One surprising feature of the theory is that 
interpolation pairs occur not infrequently. 

In general, interpolation unitaries can be hard to work with. However, in the 
special case where tp and r] arc s-elementary wavelets (also called MSF-wavelets 
with phase 0, or wavelet-set wavelets), the interpolation unitary takes the form of 
a composition operator with measure-preserving symbol. Every pair of wavelet sets 
gives rise to a measure preserving transformation on the underlying measure space 
in a natural way, called the interpolation map determined by the pair. The interpo- 
lation theory for such wavelets and their associated wavelet sets is, in many concrete 
cases, computable by hands-on experimental paper-and-pencil computations. This 
permits experimentation in the form of testing of hypotheses in potential theorems 
for more general types of wavelets. The simplest case is where a pair of wavelet 
sets has the property that the measure preserving transformation a is an involution 
(i.e. a o a = id). In this case the composition unitary has square /, so the pair of 
wavelets is indeed an interpolation pair. The pair of wavelets sets is, by analogy, 
called an interpolation pair of wavelet sets. More generally, an interpolation family 
of wavelets (and analogously, of wavelet sets) is a finite (or even infinite) family 
of wavelets for which the associated family of interpolation unitaries (interpolation 
maps) forms a group. 

It is appropriate to give a bit of background and history that serves to indicate 
why interpolation pairs of wavelet sets are relevant to the theory of wavelets. More 
exposition on this, including specific details and statements of theorems involved, 
can be found in the semi-expository articles [La2] , [La3] , [La4] , and [La5] . In [DL] , 
for any interpolation pair of wavelet sets (E, F), the authors constructed a 2 x 2 
complex matrix valued function (called the Coefficient Criterion, see [DL], Proposi- 
tion 5.4, and also [La2], section 5.22) which specifies precisely when a function / on 
R with frequency support contained in i?UF is an orthonormal wavelet. If supp{f} 
is contained in E L) F, then this criterion shows that / is an orthonormal wavelet 
iff this matrix valued function is a unitary matrix (a.e.), and it is a Riesz wavelet 
iff it is an invertible matrix (a.e). Moreover, it shows that a Parseval frame wavelet 
(resp. Riesz frame wavelet) with frequency support contained in £?UF is necessarily 
an orthonormal wavlet (resp. Riesz wavelet). It follows that the set of orthonormal 
wavelets with frequency support contained in the union of £' U F, where {E, F) is 
an interpolation pair of wavelet sets, is pathwise connected in £^(M). The set of 
Fourier Transforms of this set is also connected in the C°° norm on the frequency 
space. 

These results were the main motivating factor in posing the first open problem 
discussed in [DL] , namely the question of whether the set of all orthonormal wavelets 
in £^(M) is norm-pathwise connected. This was the same problem that was posed 
completely independently by G. Weiss and his research group in [HWW1],[HWW2] 
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for different reasons. Their reasons included the interesting discovery that certain 
wavelet sets (rather, the associated MSF wavelets) could be "smoothed" in a con- 
tinuous fashion to obtain wavelets that were continuous in the frequency domain. 
It turned out that our operator-interpolation approach, in certain key cases, was 
equivalent to the "smoothing" approach of G. Weiss, and the cases involved in- 
cluded the derivation of Y. Meyer's classic family of wavelets that are compactly 
supported and continuous in the frequency domain. Exploring common interest in 
the relationships between smoothing of a wavelet set on the one hand, and operator- 
theoretic interpolation between a pair of wavelet sets on the other hand, and the 
general "connectedness" problem that was motivated by both approaches indepen- 
dently (as described above) let to the formation of the WUTAM Consortium (short 
for Washington University and Texas A&M University) and the joint work [Wut] of 
the consortium, in which the connectedness problem was shown to have a positive 
answer for the case of MRA wavelets. 

The basic idea behind operator-interpolation is elementary. If x and y are ele- 
ments of a vector space V , we say that a vector z is linearly interpolated from x and 
y if 2; is a convex combination of a;, y. More generally, it is convenient to allow arbi- 
trary linear combinations. So the set of vectors interpolated by x and y is the linear 
span of X and y. More generally, we can say that z is linearly interpolated from a 
collection !F of vectors if z is a linear combination of vectors from the family. And 
more generally yet, if the vector space F is a left module over some operator algebra 
T> we can consider linear combinations from with coefficients that are operators 
from X>, called modular linear combinations. If 2; is a modular linear combination of 
then we say that z is derived from x and y by operator-theoretic- interpolation 
(or operator- interpolation for short). In the case of wavelets, the operator algebra T> 
is the von Neumann algebra of all bounded linear operators acting on C'^iW'') that 
commute with the dilation and translation unitary operators for the wavelet system. 
When we conjugate this with the Fourier transform (which is unitary), so we are 
working in C^iW^) as the frequency space, and if we denote this conjugated algebra 
by P, then T> is an algebra of multiplication operators on £^(M"). In particular, it is 
a commutative algebra. If a wavelet 77 is a modular linear combination of wavelets 
ip and V with coefficients which are operators in P, then we say r] is derived by 
operator-interpolation between ij) and u. Not all modular linear combinations of 
ip and V are orthonormal wavelets. They are all Bessel wavelets to be sure, but a 
certain unitarity condition needs to be satisfied to be an orthonormal wavelet. Let 
V^^ be the interpolation unitary from ip to v. If A and B are operators in D and 
rj = Atp + Bu, the necessary and sufficient condition for 77 to be an orthonormal 
wavlet is that the operator U := A + BV^^ needs to be unitary. (More generally, for 
a frame wavelet the criterion is that U must be surjective, and for a Riesz wavelet 
U must be invertible.) 



4 



XIAOFEI ZHANG AND DAVID R. LARSON 



The reason that interpolation pairs of orthonormal wavelets are special is that 
if v) is an interpolation pair, and if A and B are operators in D such that 
A * A + B * B = I, then under certain circumstances Aip + Bv can also be an 
orthonormal wavelet. In particular, li 9 G [0, 27r] is arbitrary, then r/ := cosOip + 
isinOu is an orthonormal wavelet. Indeed, in this case the operator U above is just 
cosOI + isinOV^, and since {V^'^)'^ = I it follows that UU* = I, so U is unitary, as 
required by the criterion. Letting 9 vary continuously it follows, in particular, that 
i/j and u are pathwise connected via a path of orthogonal wavelets. 

For wavelet set wavelets, i.e. MSF wavelets with phase 0, more is true, and the 
operator-algebraic geometry involved is fairly rich. For any pair of wavelets sets E 
and F, with associated MSF wavelets ipE and ijjp, the interpolation unitary V^^ 
normalizes the von Neumann algebra D in the sense that V^^^(V^)* = V. This 
was proven in Chapt 5 of [DL], and was a key result of that memoir. (We note 
that it is an open question (see [La2], Problem 4) as to whether arbitrary interpo- 
lation operators (i.e. for non-wavelet set wavelets) normalize V.) The reason that 
interpolation pairs of wavelet sets, are even more special than general interpolations 
of wavelets is the following: Firstly, {E, F) is an interpolation pair of wavelet sets 
if and only if the pair of wavelets (V'EjV'f) is an interpolation pair of wavelets. 
And secondly, since in this case the interpolation unitary normalizes V , and since 
(y^^'^-p')2 = /, it follows that the set of operators {A + BV^^ \A,B €V} is closed 
under multiplication and is in fact a von Neumann algebra. (For a more general 
interpolation pair of wavelets whose interpolation operator normalizes V the same 
thing is true.) Since the unitary group of a von Neumann algebra is pathwise con- 
nected in the operator norm, the interpolated family of wavelets is also connected. 
Much more is true. Since the elements of V are multiplication operators in a cer- 
tain family (the dilation-periodic operators), if we write A = Mf and B = Mg 
we obtain {A + BV^^)iI)e = f'lpE + fi'V'-F- Thus the outcome is an actual formula 
as well as a criterion (the Coefficient Criterion mentioned above), for constructing 
all orthonormal wavelets whose frequency support is contained in the union E U F 
of a given interpolation pair (E, F) of wavelet sets. By choosing / and g appro- 
priately, so they are continuous and vanish on the boundary of the union E U F, 
and agree on E r\ F, and satisfy the unitarity condition referred to above, one can 
obtain wavelets in this fashion which are smooth in the frequency-domain. Not 
all interpolation pairs (E, F) can be so smoothed. But some can. As alluded to 
above, Y. Meyer's famous class of orthonormal wavelets which are continuous and 
compactly supported in the frequency domain can be derived in this way from a 
special interpolation pair of wavelet sets: namely the pair E = [— x)'-'[t"'T") 
and F = [— ^)U[^,^). Even for cases in which smoothing cannot work (for 
instance, ELiF may have too many boundary points), the operator algebra involved 
can be interesting. 
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The purpose of this paper is to provide solutions to two related problems con- 
cerning dyadic orthonormal wavelet sets in the line. One problem asked whether 
a certain containment relation for an interpolation map implies that the associated 
pair of wavelet sets is an interpolation pair. This was posed by the second author in 
a VIGRE seminar course at Texas A&M several years ago. We answer this question 
affirmatively, and we also observe that the analogous result does not hold for a gen- 
eral interpolation family of wavelet sets. The second problem was posed by D. Han, 
who asked whether the equality of the congruence domains of an interpolation map 
and its inverse implies that the associated pair of wavelet sets is an interpolation 
pair. We were able to give a counterexample to this problem. Our work on this 
interesting problem motivated a useful lemma on constructing wavelet sets, which 
is apparently different from those methods which have appeared in the literature 
to date, and which is used in this counterexample as well as in the construction of 
several wavelet sets concerning some related questions, and also some wavelet sets in 
the plane. Much of the work presented in this article is material from the doctoral 
dissertation of the first author [Zh], which has not appeared elsewhere. 

Much work has been accomplished on the topic of wavelet sets since the mid 
1990's. We have outlined some of the background and history in the opening para- 
graph of Section 5, for the interested reader. Our main results in this paper are for 
the special case n = 1 with the dilation scale factor 2 and integer translates (the 
dyadic case). In Section 5, we apply one of our techniques to the construction of cer- 
tain dyadic wavelet sets in the plane. Based on our results, some further directions 
are suggested for higher dimensions. 

2. Two Problems 

A dyadic orthonormxil wavelet is a function ip G (Lebesgue mesasure) with 

the property that the set { 2 2 ^(2" • —I) \ n,l € Z } forms an orthonormal basis for 
C'^ (M) . More generally, if A is any real invertible n x n matrix, then a single function 
il) G L^(M") is an orthonormal wavelet for A if 

{|deU|i^(A" • -Z) I n G Z, Z G Z^")} 

is an orthonormal basis of L^(M"). If A is expansive (equivalently, all eigenvalues of 
A are required to have absolute value strictly greater than 1) then it was shown in 
[DLSl] that orthonormal wavelets for A always exist. 

By the support of a measurable function we mean the set of points in its domain 
at which it does not vanish. By the support of an element / of jC^(M), we mean the 
support of any measurable representative of /, which is well-defined in the measure 
algebra of equivalence classes of sets modulo null sets. By the frequency support of 
a function we mean the support of its Fourier Transform. 
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Let T denote the n-dimensional Fourier transform on £^(M") defined by 



(^/)(^) ■■ 



(27r)2 h 




for all / € L^(]R"). Here, sot denotes the real inner product. A measurable set 
E C is a wavelet set for A if 



is an orthonormal wavelet for A. 

A sequence of measurable sets {En} is called a measurable partition of E if 
E = EA{\^^En) is a null set and En fl Em has measure zero if n 7^ m, where 
A denotes the symmetric difference of sets. Measurable subsets E and i*" of M are 
called 277 -translation congruent to each other, denoted by E ~27r if there exists 
a measurable partition {En} of E, such that {En + 2mr} is a measurable paritition 
of F. Similarly, E and F are called 2-dilation congruent to each other, denoted 
by E2^ F, if there is a measurable partition {En} of E, such that {2'^En} is a 
measurable partition of F. A measurable set E is called a 2Tr -translation generator 
of a measurable partition of R if {E + 2mr}neZ forms a measurable partition of 
M. Similarly, a measurable set F is called a 2-dilation generator of a measurable 
partition of M if {2"F}„g2 forms a measurable partition of R. 

Lemma 4.3 in [DL] gives the following characterization of wavelet sets, which was 
also obtained independently in [FW] using different techniques. 

Let E CM. be a measurable set. Then E is a wavelet set if and only if E is both 
a 2tt -translation generator of a measurable partition o/M and a 2-dilation generator 
of a measurable partition o/M. 

Again from [DL], suppose that E, F arc wavelet sets, and let a : £' ^ i*" be the 
bijective map (modulo null sets) implementing the 27r-translation congruence. Then 
a can be extended to a bijective (modulo null sets) measurable map on M by defining 
(t(0) = and a{s) = 2"cr(2-"s) for each s G 2"E, n G Z. This map is denoted by 
(7^ and called the interpolation map for the ordered pair of wavelet sets {E, F). 

An ordered pair of wavelet sets {E, F) is called an interpolation pair if o'^ocr^ := 
(o"^)^ = id^. In general, an interpolation family of wavelet sets is a family of 
wavelet sets such that { | F G .F } is a group under composition of maps for some 



Question A. Let E, F be wavelet sets. Does CF^i^ U F) C E U F imply that 



E 



{a^f = idm? 



We answer Question A affirmatively. Wc then give an elementary example that 
shows that this result need not hold for a triple of wavelet sets; However, there is a 
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natural modification which does make sense for n-tuples, and we pose it as an open 
question (Question C). 

Given a pair of wavelet sets, it is not obvious at all upon initial inspection whether 
they actually form an interpolation pair. And constructing interpolation pairs can 
be hard. A basic problem from [DL] that still remains open is the question: Given 
an arbitrary dyadic wavelet set E in M}, is there necessarily a second distinct wavelet 
set F such that {E, F) is an interpolation pair? Partly to address this problem, and 
partly for intrinsic interest, Han introduced and studied properties of congruence 
domains in [Hani] and [Han2]. Given a pair of wavelet sets {E,F), the domain 
of 2TT -congruence of cr^, denoted by V^, is the set of all points s G M such that 
(T^(s) — s is an integral multiple of 2tt. There is a close relation between this and 
the interpolation map of a pair of wavelet sets. Han asked the following question. 

Question B. Let E,F be wavelet sets. Does = Vp imply that (iTg)^ = id^? 

If {E,F) is an interpolation pair, then it is easily verified that the domains of 
27r-congruence of {E, F) and {F, E) are the same. So the above question just asks 
if the converse is true. In many cases it is true. But it is not universally true. 
We answer Question B negatively by constructing a counterexample. The key is a 
special lemma on constructing wavelet sets, which is also used to build examples for 
several other questions. 

3. SOLUTION TO QUESTION A 
The following theorem provides an answer to Question A. 

Theorem 1. Let {E,F) be wavelet sets. The following two statements are equiva- 
lent: 

(i) {E,F) is an interpolation pair, 

(ii) a^{E UF)CEUF. 

Proof, (i) =^ (ii). Suppose that {E,F) is an interpolation pair. Then (cr^)^ = 
id^. Since {o-p)~^ = ap, it follows that = af. Observe that o-p{E) = F, 
af{F\E) = E\F, and so 

a^{E UF) = al{E) U (tI{F\E) = F U (jf{F\E) = F U E\F = E U F. 

(ii) ^ (i). Suppose that a^{E U F) C E U F. Since a^{E) = F and erf is 
bijective (modulo null sets), we must have ap{F\E) C E\F. We will prove this by 
way of contradiction, and a diagram is included at the end of the proof to help in 
navigating between E and F, with arrows representing a p. 

Assume that {E,F) is not an interpolation pair. Let Gq = {s G E \ (erf )^(s) ^ 
s}. Then Go C E\F is Lebesgue measurable and has positive measure. Since 
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{Gq n {F — 2mT) }nez forms a measurable partition of Gq, it follows that Gq fl 
{F — 2ni7r) has positive measure for some ni G Z. Denote this set by Gi. Then 

'^s(^i) = Gi + 2ni7r C F\E. Since is 2-dilation congruent to F\E, following 
the similar discussion, there exists a measurable subset G2 of Gi with positive 
measure, such that 2~''i(G2 + 2ni7r) C E\F for some integer ki. Then, there exists a 
measurable subset G3 of G2 with positive measure such that a^{2~^^ {Gs + 2ni7r)) = 
2~''^{Gs + 2ni7r) + 2n27r for some integer n2- Thus, 

{'^E?{G3) = cj^iGs + 2ni7r) = 2^^ ■ a^i2-^' (G3 + 2ni^)) 

= 2^=1 • (2-^=1 (G3 + 2rai7r) + 2n27r) = G3 + 2rai7r + 2^=1 • 2n2Tr C 

Since both G3, G3 + 2ni7r + 2^^ • 2n2'7r C E\F and they are distinct by assumption, 
we must have rii + 2^'^ ■n2 ^ Ij, which implies that ki < —1. Similarly, there exists a 
measurable subset G4 of G3 with positive measure, such that a^{G4 + 2ni7r + 2^^ ■ 
2n27r) = G4 + 2ni7r + 2''i • 2n27r + 2n37r for some integer n^. Then, 

(a^)2(2-^'i(G4 + 2ni7r)) = a^{2-^^{Gi + 2ni^) + 2n27r) 
= 2-''i • C7^(G4 + 2ni7r + 2^^ • 2n27r) = 2-'''i • (G4 + 2nnr + 2*^1 • 2ra27r + 2n37r) 
= 2"*^! (G4 + 2ni7r) + 2n27r + 2'''^ • 2n37r C E\F. 

Since both 2~''^ (G4 + 2ni7r) and 2~''^ (G4 + 2ni7r) + 2n27r + 2~''^ ■ 2n37r are contained 

in E\F and n2 + 2"^^ • n3 must be an integer, we must have n2 + 2~^^ • 77.3 = 0. 
Then, since maps both G4 and G4 + 2ni7r + 2^'^ ■ 2n2'^ to G4 + 2ni7r, G4 and 
G4 + 2ni7r + 2^'^ ■ 2n2vr must be the same, then (o"^)^(s) = s,Vs € G4 C Gq, 
contradicting the assumption. Thus, (cg)^ = id^. almost everywhere. 

E\F : Go 2-'=i(G2 + 2ni7r) 

I I 

F\E : Gi + 2ni7r 2-'=i (G3 + 2nnr) + 2n2n 

I I 

: G3 + 2ni7r + 2*^1 • 2n27r 2'''^ (G4 + 2ni7r) + 2n27r + 2'''^ ■ 2n37r 

I I 

: G4 + 2ni7r + 2*^1 • 2n27r + 2n37r 

□ 



Let JT be a general interpolation family of wavelet sets, and fix Ei € Then 
for arbitrary £"2 , -Es G we have for some iiJ4 € Thus 

(£;3) = oaf = = £?4. it fillows that (Uce^^) ^ Uce^^ 

for each Ex G J'^. However, the converse may be false as the following example shows. 

Example 2. Let (E, F) be an interpolation pair of wavelet sets and suppose 
G C ii^ U F is a wavelet set contained in the union which is distinct from E and 
F. Observe that E\F is both 27r-translation congruent and 2-dilation congruent to 
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F\E, thus G C E Li F must contain E r\ F to he a wavelet set. Since = id^ 
on E Ci G and cr^ = on E\G, ((Jg)^ = id^. Similarly, (cTp)^ = ids.. Then, 
° ~ '^F 7^ '^d.-R, which implies that {id^,a^,a'^} is not a group under com- 
position of maps. Hence {E, F, G} is not an interpolation family. However, no- 
tice that {E\G, G\F, EOF, G\E, F\G} forms a measurable partition of U F. 
E\G is both 27r-translation and 2-dilation congruent to G\E and F\G is both 27r- 
translation and 2-dilation congruent to G\F which is contained in £' fl G. We have 
a^{E UFUG) = a^{E UF) = EUF, and 

a^iEUFuG) = a'i{EUF) = a'iiE)Ua'iiF\E) 

= GU (7%{G\E) U (T%{F\G) = G U E\G U F\G = EUF. 

□ 

Observe that in the above example, if let H = E\G D {E F) U F\G, then 
H Q EUF is also a wavelet set and cr^ = o cr^ = o cr^. Then {E, F, G, H} is 
an interpolation family since {a^ = idu, o"^, cr^, cr^ } is actually isomorphic to the 
Klein four group. This motivates the following question which seems difficult and 
which we pose as an open problem: 

Question C. Let ^ be a finite collection of wavelet sets. Fix Ei G J^. Sup- 
pose cr^^ (Ufie.^ ^) contained in IJ^e.^ ^ E\ & J^. Is .7^ a subset of a 
finite interpolation family of wavelet sets? That is, can !F always be extended to a 
finite interpolation family? 

4. A COUNTEREXAMPLE TO QUESTION B 

Before giving a counterexample to Question B, we present a characterization of 
interpolation pairs of wavelet sets and equality of 27r-congruence domains of asso- 
ciated interpolation maps. This is actually the motivation that led us to find a 
counterexample . 

Given a pair of wavelet sets (^E^F^y there exists a measurable partition {-E'nlnGZ 
of E such that {£"„ + 2mT}nei is a measurable partition of F. Similarly, there exists 
another measurable partition {E^}keZ of E such that {2^E''}keZ is another mea- 
surable partition of F. Denote En H E^ by En^k, then {En,k}n,kez is a measurable 
partition of E, and both {En^^ + 2n7r}„_feg2 and {2^ En^k}n,k&. are measurable par- 
titions of F. Observe that if E^^k has positive measure, then n = whenever A; = 0, 
and En^k H iP^^Fm^i — 2mr),n, k,m,l G Z\0 forms a partition of E\F. 

Theorem 3. Let E, F be wavelet sets, and let {-E'n,fe}n,feez be as defined above. 
Then: 



10 



XIAOFEI ZHANG AND DAVID R. LARSON 



(i) {E, F) is an interpolation pair if and only if E^^k ^ (2'-£^m,/ ~ 2n7r) with positive 
measure for some nonzero integers n, k, m, I implies that n+2' -m = 0. Furthermore, 

for such a set, we also have k = —I. 

(a) = Dp if and only if En,k H {2^Em.i — 2mr) with positive measure for some 
nonzero integers n,k,m,l implies that [n] + 1 = [m]. Here, [n] denotes the smallest 
integer k such that 2*^ • n G Z. 

Proof, (i) Suppose that {E, F) is an interpolation pair. Then (o"^)^ = id. Since 



{a^)\E^^kni2^Em,i-2mr)) = a^{{En,k + 2mr) n2^Em,i) 

= 2^ ■a^{2-\Er,,k + 2n7r)nE^^i) 

= 2^ ■{{2-^{En,k + 2mT)nEm,i) + 2rmT) 

= + 2n7r + 2' • 2m7r) n 2' + 2m7r) , 



it follows that if En^k H {2 E^^i — 2mr) has positive measure, then 2mr + 2 • 2rmr = 0, 
i.e., n + 2' • m = 0. 

Conversely, suppose that En^k H {2^Em,i — 2n7T),n,k,m,l ^ having positive 
measure implies that n + 2' • m = 0. Observe that {En,k H (2'-^^^/ — 2n'K)}n,k,m,l^o 
forms a measurable partition of E\F. Then, a similar argument shows that ((Jg)^ = 

Furthermore, if E^^k ^ (2'-E'm„/ ~ 2mr) having positive measure implies that n + 
2' • m = 0, then {Em,i + 2m7r) n 2^^En^k = 2"' • {En^k n {2'-E„^ - 2n7r) ) also has 
positive measure. Since Ejn,i + 2m7r C F and 2~^En^k £ F only if = —I, we must 
have k = —I. 

(ii) Observe that 



= U i[J^'E^,k) U (j2^{EnF) 

n,k^O j>[n] j 

= U ( U 2^{En,k n 2-\Em + 2m7r))) U \j2^iE n F), 

n,k,mj^O j>[n] j 
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and 

^1 = U ( U 2'(^n,fe + 2n7r)) U [j2\E n F) 

n,fey^O j>[n] j 

U ( U + 2n7r) n 2'E^,i)) U (J2^'(^ n F) 

n,k,m,l^O j>[n] j 

= U ( U 2''^'(^m,; n 2-'(i?„,fc + 2n7r))) U U2'(^ n F) 

n,k,m,l^O j>\n\ j 

= (j ((j 2^'+'(^n,fc n 2-\E,n,i + 2rmr))) U \J2^ (E D F) 

n,k,m,ljtO j>[m] j 

= U ( U 2'^'{En,kn2-''iEm + 2m7T))) U \j2^{EnF), 

where U denotes disjoint union. The next to last equahty comes from changing 
indices. Thus, = Vf imphes that [m] + k = [n], if E^^k n 2^'^{Em + 2m7r) has 
positive measure, or equivalently, if E^^i fl (2'^£'„ ^ — 2m7r) has positive measure for 

some I / 0, since 2^{En.k n 2-^{Em + 2m7r)) - 2m7r = (2'=£;„,fc - 2m7r) ^ Em = 
{2^En,k — 2m7r) fl (\JiEjn,i) ■ The converse direction can be shown by reversing the 
above discussion. □ 



The basic idea we use is to construct two certain measurable subsets of M, which 
are both 27r-translation and 2-dilation congruent to each other, yet the interpolation 
map restricted to the union of 2-dilates of which does not have the property that 
its square equals the identity map, and then to construct the remaining pieces of 
wavelet sets for these two sets so that the congruency domains match up. This type 
of approach was used by D. Speegle in [SI] and Q. Gu in [Gul], and a necessary 
and sufficient condition has also been given for a measurable set being contained in 
some wavelet set in [IP]. However, the methods and the constructions we use in the 
present paper are completely different from those in [SI], [Gul] and [IP]. 



Theorem 4. The answer to Question B is no. 



Before proving Theorem 4, we require a technical lemma. This will also be used 
in constructing several other examples in the next section. 

Lemma 5. Let E,F C. R" be Lebesgue measurable sets with finite positive mea- 
sure. If there exist n\,n2,k\,k2 & such that modulo null sets, 

2^^ F (Z E + 2niTT and F + 2n27r C 2^^F, 

and {E + 2n27r) n 2^^F is a null set, then there exists a Lebesgue measurable set 
G C. [E + 2n27r) U 2^iF such that G is both 2n -translation congruent to E and 
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2-dilation congruent to F. In fact G can be taken as: 

oo 

G = (J 5^(2^=1 F\2'=i-'=2(E + 2n27r)) 

i=0 

oo 

U (S + 2n27r)\(lJ 2^=2-^1 . S'{2''^ F\2''^~''^{E + 2n2Tr))), 

i=l 

where S{x) = 2^^-^''{x + 2(n2 - ni)7r), Vx G M". 
Proof. By hypothesis, 

2fei-fc2(£; + 2n27r) C 2*^iF C £; + 2ni7r. 

Since ^i{2''^-^^{E + 2n27r)) = 2^^-^^ ii{E) < fi{E + 2ni7r) = //(S), h - k2 < 0. If 
fei = A:2, then + 2n27r C 2''^F, contradicts the hypothesis that (E + 2n2vr) n 2''^F 
is a null set. Thus, fci — ^2 < 0. Construct a sequence of measurable sets 
as follows. Let 

Go = 2''^F\2''^-''^{E + 2n2Tr) 

= 2^^F \ 2^^-^^{E + 2ni7r + 2(n2 - ni)7r) <^ {E + 2nnr) \ S{E + 2nnr). 

Let Gi = 5* (Go) C S'{E + 2nnr)\S'+^{E + 2ni7r), for each i € N. Notice that 
the measure of Gj is bounded by 2*'('^i~'^2) . yrt(£;), which will approach as i ap- 
proaches infinity. Furthermore, S{E + 2ni7r) C E + 2ni7r, and it follows that the 
Gi's are disjoint and [JZoGi C 2''^F. Also by definition, W^^^Gi C S{E + 2ni7r) 
and 2'=2-'=iGi+i = Gj + 2(n2 - ni)7r, > 0. Let 

CXD OO 

G = {\jGi) U (£; + 2n27r)\(|j2'=2-*^iG,). 

1=0 i=l 

Then, since 2'=i~^2((s + 2n27r)\(U=i 2'^2-fciGr.)) = + 2ni^)\(Ul^i G,), and 
{Gi}°Zo U {'S'(-E + 2ni7r)\(Uj^i Gi)} constitutes a measurable partition of 2*^1^, it's 
clear that G is 2-dilation congruent to F. On the other hand, since 

oo oo 

(J 2'=2-fciGi = U(Gi + 2(«2 - ni)7r), 

i=l i=0 

G is 27r-translation congruent to □ 

Proof of Theorem I Let Ei = [f|7r, f|7r), ^2 = [f vr, f 7r)+127r, = [f tt, f 7r)+87r 
and ^4 = [f^,f^) + 96vr. Let Fi = [f7r,f^),F2 = [Ivr, f|7r) + 67r, F3 = 
[^TT, ^tt) + 167r and F4 = [^vr, ^vr) + 24tt. Notice that E^i , ^2 , -^3 , ^^4 are both 
27r-translation congruent and 2-dilation congruent to disjoint pieces of [— 27r, — tt) U 
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E1UE2UE3UE4 F1UF2UF3UF4 ~2,r [^,7r), 

El U E2 U E3 Li E4^ Fi U F2 U F3 U F4 2- 

. TT TT 7 TT 7 TT 

f^+128'^+16^^[r+128'r+64)- 
Let G be the measurable set determined by Lemma 5 and the fohowing two con- 
tainment relations: 

1 /r 71" N r TT 7 TT , .7 Vr „ r„ TT , 

[0,i) + 8, c 8.(|,,, + :^)U|, + i;;.I.+ ^)U[Z, + iL,2.)). 

Then, G is both 27r-translation congruent to [0, j^) and 2-dilation congruent to 
[tt, tt + ^) U [tt + ^, Itt + ^) U [|7r + ^, 27r). Hence, 

:= [-27r, -tt) U G U U £;2 U £^3 U £^4, 
F := [-27r, -tt) U G U Fi U £2 U F3 U F4 
are both wavelet sets. The associated interpolation map cr^ is defined by: 

s + 67r if s G Fi 



s + 127r if s € F2 
s - Stt if s G F3 



SOtt if s G F4 

if s G [-27r, -tt) U G. 



Straightforward computation shows that 

Vl == (-00, 0) U (J 2'^G U (J (Fi U Jf2 U ^F3 U ^F4) 



A;>0 



and 



Vf == (-00, 0) U (J 2^'G U U (^2 U ^F4 U jFi U ^F3). 



fe>0 



Observe that F2 = ^F2, 5F4 = |F4, |Fi = Fi, and IF3 = ^E^, which implies that 
= Vp. However, since 

(a^fiEi) = a^{Ei + Gtt) = a^{^E2) = ^a|(F2) = ^(Fs + 127r) = Fi + 127r, 
(F, F) is not an interpolation pair. □ 



It turns out that if one of F, F is the Shannon set, then the equality of integral 
domains does imply that (F, F) is an interpolation pair. 



Proposition 6. Let F, F be wavelet sets. If E = [-27r, -7r)U[7r, 2tt), then P| = Vf 
implies that (F, F) is an interpolation pair. 
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Proof. Let = E H {F — 2mr), n G Z. Then {£?n}nez is a measurable partition of 
E and {En + 2mr} is a measurable partition of F. Let E~ = EnC) [— 27r, —Tr),E^ = 
Eji n [tt, 27r). Then, we have the following diagram. 

E : 



— 277 — TT TT 27r 



jp . ■■■ EZi-2'K Et2-4=T^ Eq Eti-2'K E^+27r E^ E^^tt E^+2Tr ■ ■ 

A A \ \ \ \ A A A 

— 47r —Sir —2-K — tt tt 27r Stt 47r 



Observe that 

Vl= U 2i{E-UEt) U [j2^E,. 

n7^0j>[n] j& 

Then 

n[0,7r) = U 2^Et U \J2^E+ 

n^O,[n]<j<0 j<0 
= U '^'^t U ^^^2 U ^^^4 U \Et^ U i^+e U \e+, U ii^+s U . • ■ , 

here E±n ■= E^ U E^n- ^| n [tt, 27r) = [tt, 27r). Similarly, 

Pf = U 2^((^-U^+) + 2n7r) U U2''^o. 

Then, following from the diagram, 

n [0, tt) = (^f + 27r) U U^^-^o"' 

P| n [tt, 2tt) = E+{J {\J2^{E^ + 2tt) n [tt, 27r)) U ^(£;2- + 47r) U + Stt) U • • • . 

Comparing P| n [0, vr) and n [0, vr), we have 

(1) Er+27r = ^£;i2u^^i4uj^i4u^£^i6u^£^i8u^^i8u---. 

Then, it is easy to see that (Uj>i 2-'(-Ej~ + 27r)) fl [tt, 2it) = Une2Z\o ^n- Comparing 
n [tt, 27r) and n [tt, 27r), we have 

(2) U ^ ^iE-+i7r)iJ^iE-+87r)ij--- . 

ne2Z+l 

By symmetry, we also have 

(3) £;+ -27r = \Eg,ij^Eg^U^Eg^U^E^,U^E^,U^E^,iJ ■■■ , 

(4) (j E- C ^{E+,-A7r)U^{E+,-87r)U--- . 

ne2Z+l 
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(1) and (4) imply that + 27r = ^E^2j ^^'^ measure zero for each odd 
integer n except 1, and has measure zero for each even integer n except —2, 0. 

(2) and (3) imply that E^-^^ — 2ti = ^-E^, and E'^ has measure zero for each odd 
integer n except —1, and E^^ has measure zero for each even integer n except 0,2. 
Thus, -F must have the following form 

F = {E+2 - 47r) \JE^{J {Eti - 27r) U {E^ + 27r) \J E+ \J {E^ + 47r). 

Using the fact that E^ +2it = ^E^2 ^^'^ ^-i — 27r = ^-E^, it is easy to verify that 



5. Some examples of wavelet sets 

Wavelet sets are useful as examples and counterexamples. Many examples of 
them in the real line were given in [DL] exactly for that purpose, for experimenta- 
tion and testing hypotheses, and many open questions still remain in that setting. 
The existence of wavelet sets in the plane, and more generally in M", was first proven 
by Dai, Larson and Speegle [DLSl] in the summer of 1994 during a course at Texas 
A&M taught by the second author using the manuscript of [DL] as a text. The 
proof in [DLSl] was abstract and covered dual congruence results for more general 
types of dual-dynamical systems. But given the definitions, the proof was construc- 
tive, implicitly giving an algorithm (a casting- outward technique) for constructing 
examples of wavelet sets for arbitrary expansive matrices and full rank translation 
lattices on M". However, this method produced only wavelet sets that were un- 
bounded and had as a limit point, and were not very well described by diagrams 
or pictures. In other words, they were wavelet sets but not "nice" wavelet sets. In 
1995 concrete examples of dyadic wavelet sets in the plane, which were bounded and 
bounded away from 0, and could be nicely diagrammed, were constructed by Soardi 
and Weiland [SW]. At about the same time, Dai and Larson constructed two 
dyadic wavelet sets (denoted the "four-corners set" and the "wedding cake set"), for 
inclusion in a final version of [DL] in response to a referee's suggestion. In 1996-97, 
a number of authors constructed many other concrete examples of wavelet sets in 
the plane, and in higher dimensions. Included are the wavelet sets computed and di- 
agramed in the articles [BMM], [BLl], [BL2], [DLS2], [Za], [C], [GH], [Gu2], [Han2], 
[S2]. The last three are the Ph.D. theses of Gu, Han and Speegle, respectively. 
Importantly, Baggett-Medina-Merrill [BMM], and Benedetto-Leon [BL2], indepen- 
dently found two interesting completely different constructive characterizations of 
all wavelet sets for expansive matrices in M". Open questions remain, especially 
questions concerning the existence of wavelet sets with special properties, and algo- 
rithms for constructing special classes of such sets. 

In this section, we will apply the techniques introduced in Lemma 5 to construct 
an unbounded symmetric wavelet set and counterexamples to two related questions 
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on wavelet sets, in addition to new constructions of some known wavelet sets. 
Example 7. Let Gn = [^^^tt, ^^t+T^ 7r) + (2"+^ - 2)tt, n>0, then 

_ 1 on+l _ 1 

G„ C [(2«+2 _ 2)n, (2"+2 _ 1)^) = 2«+i [vr + ^^tt, vr + ^^tt). 

Observe that U„>o(G'n " (2"+' - 2)7r) = [f ,7r) and Un>0 2^<^n C Un>o[^ + 
^vr, vr + ^5i?i7r) = [fyr, 27r). Let Eo = [0, f ), let = [tt, 27r) \ {(j'^^^^Gn). 
Notice that since Fq D [tt, |7r), 

£;o + 27r C 2Fo and ^Fq C £;o, 

and (£'o + 27r) n = 0. Let Go be the set obtained by applying Lemma 5. 
Then ij^o^" both 27r-translation congruent to [0, tt) and 2-dilation congruent to 
[tt, 27r). Hence, by symmetry, — (Un=o^") ^ (Un=o^") unbounded symmetric 
wavelet set. □ 



We note that the first example of an unbounded symmetric wavelet set was given 
in Proposition 2.14 of [FW]. The construction involved some significant computa- 
tions and explanation. Example 4.5 (xi) of [DL] gave a different unbounded wavelet 
set whose construction required little explanation, but it was not symmetric. 

The following is a counterexample to a question posed by the second author in 
[Lai]: 

Let E be a wavelet set, and suppose that G C 2E U E U ^E is also a wavelet set. 
Is {E, G) an interpolation pair? 

Example 8. Let E = [-2-k, -tt) U [vr, 2-k), and G C 2E U E U ^E = [-An, -f )U 
[^, 47r) be a wavelet set. Suppose that Gi = GO [3tt, Att) has positive measure, then 
Gi - 27r C [7r,27r) C E. Vs G Gi - 27r, (ag)2(s) = (s + 27r) = 2 • a'^{^s + tt) = 
2 ■ (is + TT + 2kTr) = s + 27r + Akir, for some k e Z. Thus, (erf 7^ s for each 
s G Gi — 27r. Therefore, {E, G) is not an interpolation pair if Gfl [Stt, 47r) has positive 
measure. Since 

^[7r,^7r) C [0,7r) and [0, tt) + 27r C 2[7r, ^tt), 

and i[7r, |7r) n ([0, vr) + 2tt) = 0, by Lemma 5, there exists Gq C [^, |7r) U [27r, Svr), 
such that Go is 27r-translation congruent to [0, vr) and 2-dilation congruent to [tt, jtt). 
Let G = [-TT, -f ) U Go U [^TT, 47r). Then G is a wavelet set and G C 2E U U \E, 
but (£', G) is not an interpolation pair. □ 

The next example answers another question from [Lai]: 

Let E he a wavelet set, and suppose that G C. (E — 27r) U -E U (-E -|- 27r) is also a 
wavelet set. Is {E, G) an interpolation pair? 
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Example 9. Let E = [-27r, — tt) U [tt, 27r) and G C {E - 2Tr) U E U {E + 2tt) = 
[— 47r, — Svr) U [—2tt, 2tt) U [Svr, Att) be a wavelet set. Suppose that Gn [—tt, tt) is not a 
null set. Without loss of generality, assume that Gi = Gn[0, vr) has positive measure. 
Then, for each s G Gi - 27r C ^, (o-i)^(s) = ^^(s + 27r) = 2''' ■ af (2^(s + 27r)), for 
some k>0 such that 2'=(s+27r) € [7r,27r). Then (erf = 2-'=-(2'=(s+27r)+2n7r) = 
s + 27r + 2~'^ ■ 2mr for some n G {1, 0, -1}. Thus, (o-f )2(s) 7^ s for each s e Gi - 2it. 
Therefore, if Gfl [— 7r,7r) is not a null set, {E,G) is not an interpolation pair. Based 
on the above observation, let Go be the measurable set given by Lemma 5 and the 
following two containment relations: 

[-27r, -^tt) C [-27r, -^vr) U [-^tt, -tt), 

([_27r, -^tt) U [-^tt, -tt)) - 27r C 2[-27r, -^tt). 
Then, G := [|7r, |7r) U [tt, jTt) U [|7r, 47r) U [-f tt, -|) U Gq is a wavelet set, since 

[0,27r) = ([-27r,-^7r) + 27r)U[§7r,|7r)U([-^7r,-7r) + 27r) 



5 3 TT 7 

U [;r, -tt) U ([--tt, - J) + 27r) U ([-, 47r) - 27r), 



and 



[-27r, -tt) U [-tt, -tt) = [-27r, --^) U 2 • [--vr, --) 

^t^'^' ^^)U^-[^'r,47r)U[7r,^7r) 

is a 2-dilation generator of a measurable partition of R. However, G fl [— tt, tt) is not 
a null set, hence (cr^)^ 7^ zc/r. □ 

Lemma 5 can be useful in constructing special wavelet sets, and in particular, 
parametric families of wavelet sets. 

Example 10. Let /,m,n G N be given. Let E+{m) = [0,2-™+i7r), F+ = [7r,27r). 
Observe that 

C E+im) and F+(m) + 2^ C 2'F+, 

and 2^-F+ fl {Ej^{m) + 2'7r) = 0. Let G+{l,m) be the measurable set obtained by 
applying Lemma 5. Straightforward computation shows that 

G+(Z, m) = (^-^! TT, 2-"^+^) U [2V, 2'7r + tt). 

Similarly, let E_{m) = [-2tt + 2-"^+^, 0), F_(m) = [-27r + 2-™+!?!, -tt + 2-'"7r). 
Observe that 

F_(m) C E_{m) and F_(m) - 2"*+'^7r + 2"7r C 2"*+'^F_(m). 
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By applying Lemma 5 again, we get 

<rm+n nn 

G_(m,n) = [-2"^+"7r + 2"7r - ^^^^ _ ^ tt, -2'"+»7r + 2"7r) 

nm+n on 

U [-27r + 2-'"+V-^— —tt). 

Thus, for fixed l,m,n, G-{m,n) U G^{l,m) is both 27r-translation congruent to 
E_(m)U£;+(m) = [-2n+2-"^+'^Tr, 2""^+^) and 2-dilation congruent to F_(m)UF+, 
hence is a wavelet set. □ 



Notice that G-{1, m, n) U G+{1, m, n), Z, m, n > 1 is exactly the family of wavelet 
sets Ki^jn,n,i,m,n > 1 introduced by X. Fang and X. Wang in [FW] (see [FW] 
Example (5)). 

In the following examples, we will construct some known and new wavelet sets in 
with dilation matix A given by 2 • id]^2. It is known (see [DLSl] and [SW]) that 
£^ C is a wavelet set for 2-id^2 if and only if E is both 27r-translation congruent to 
[— TT, tt) X [— TT, tt) and 2-dilation congruent to [— 27r, 27r) x [— 27r, 27r) \ [— tt, tt) x [— tt, tt). 
Wavelet sets in higher dimensional spaces with respect to arbitrary real expansive 
matrices can also be obtained in the similar way. 

Example 11. Consider the first quadrant. Since we have the following two re- 
lations: 

[0, tt) X [0, tt) \ [0, |) X [0, |) C [0, vr) x [0, vr), 

[0, tt) X [0, tt) + (27r, 27r) C 4 • [0, tt) x [0, tt) \ [0, |) x [0, |), 

and [0, tt) x [0, tt) \ [0, f ) x [0, f ) n ([0, tt) x [0, tt) + {2tt, 2tt)) = 0. Applying Lemma 
5, we can construct a set Wi which is both 2-7r-translation congruent to [0, vr) x [0, vr) 
and 2-dilation congruent to [0, vr) x [0, vr) \ [0, ^) x [0, ^). Symmetrically, construct 
sets W2,W3,W4 in second, third and fourth quadrants, respectively. Then Wi U 
W2 U W3 U W4 is a wavelet set. Straightforward compuation shows that it is exactly 
the "four corners set" in [DLS2]. □ 

Example 12. Consider the right half plane. Since we have the following two 
relations: 

[0,7r) X [-tt,tt)\[0,^) X [-|,|) C [0,7r) x [-tt,tt), 

[0,tt) X [-7r,7r) + (27r,0) C4-[0,7r) x [-tt, tt) \ [0, |) x [-|,|), 

and [0, tt) x [-tt, tt) \ [0, f ) x [-§ , f ) n ([0, tt) x [-tt, tt)+{2tt, 0)) = 0. By Lemma 5, we 
can construct a set Wi which is both 27r-translation congruent to [0, tt) x [—tt, tt) and 
2-dilation congruent to [0, vr) x [— vr, vr) \ [0, ^) x [— ^, ^). Symmetrically, construct 
the set W2 in the left half plane. Then Wi U 11^2 is a wavelet set. Straightforward 
computation shows that it is exactly the "wedding cake set" (Example 6.6.2 in [DL] 
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and also Figure 2 in [DLS2] . ) □ 

Example 13. Consider the left-top half plane (above the line y = x in the left 
half plane). Let Ei = {{x,y)\x > — vr,y < vr, y > x}, and Fi = {(x,y)\x > — vr. y < 
TT, y > x} \ {{x, y)\x > — ^, y < ^,y > x}. Then we have the following two relations: 

Fi C El, El + (-27r, 27r) C 4 • Fi, 

and Fi fl {Ei + (— 27r, 27r)) = 0. Now use Lemma 5 to construct a set Wi which is 
both 27r-translation congruent to Ei and 2-dilation congruent to Fi . Symmetrically, 
construct a corresponding set W2 in the right-bottom half plane. Then Wi U W2 is a 
wavelet set. The diagram for this is given in Figure L □ 



r 




A 



Figure 1 Pine Tree 



Example 14. Consider the first quadrant. Let Ei = [— |, 0) x [— |, 0)U [— tt, — f tt) x 
[— TT, — Itt) and Fi = [|7r, 27r) x [|7r, 27r). Then we have the following two relations: 

F1QE1 + (27r, 27r), Ei + (47r, 47r) C 2 • Fi, 

and Fl fl {Ei + {Att, Att)) = 0. By Lemma 5, we can construct a set Wi which is both 
27r-translation congruent to Ei and 2-dilation congruent to Fi. Symmetrically, we 
can define E2,Es,E4 and ^2,^3,^4 , and construct sets W2, W3, W4 in the second, 
third and fourth quadrants, respectively. Let B = [— vr, vr) x [—tt,tt) \ ([— ^vr, ^vr) X 
[-ivr, ivr) U [fvr, vr) x [|^, vr) U [-n, - |7r) x [fvr, vr) U [fvr, vr) x [-tt, -f tt) U [[-tt, -|7r) x 
[-TT, -fvr)). Then since S U Fi U F2 U F3 U F4 = [-tt, vr) x [-tt, tt) and 25 U Fi U 



20 



XIAOFEI ZHANG AND DAVID R. LARSON 



^2 U F3 U F4 = [-27r, 27r) x [-27r, 27r) \ [-tt, tt) x [-tt, tt), VFi U 1^2 U U W4 U 5 
is a wavelet set. Computation shows that it is one of the wavelet sets introduced in 
[SW] . The diagram for this is given in Figure 2. □ 



-4pi 



-2pi -pi 




pi 2pi 



4pi 



Figure 2 A Wavelet Set in [SW] 

Remarks. (1) The idea of Lemma 5 was used in constructing a covering of R by 
symmetric wavelet sets, and this was a key to subsequent work of the first author 
with Rzcszotnik in [RZ]. 

(2) Further Directions: Interpolation maps, and interpolation pairs and more 
general interpolation families of wavelet sets, make sense and have been studied for 
matrix dilations in R". Congruence domains also make sense for matrix dilations. 
Lemma 5 in this article was stated and proved for R", and was applied to solve 
Question B in R^, and also used to study examples of dyadic (i.e. for matrix dilation 
21) wavelet sets in the plane. The results in this paper suggest some directions for 
further research. In particular, does Theorem 1 extend to matrix dilations in M", 
especially for the cases where interpolation pairs are known to exist? It might be 
useful to try to extend Lemma 5 and also Theorem 3 to matrix dilations. 
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